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Harmonic generation in the limit of ultra-steep density gradients is studied experimentally. Ob- 
servations demonstrate that while the efficient generation of high order harmonics from relativistic 
surfaces requires steep plasma density scale- lengths (L p /X < 1) the absolute efficiency of the har- 
monics declines for the steepest plasma density scale- length L p — >• 0, thus demonstrating that 
near-steplike density gradients can be achieved for interactions using high-contrast high- intensity 
laser pulses. Absolute photon yields are obtained using a calibrated detection system. The effi- 
ciency of harmonics reflected from the laser driven plasma surface via the Relativistic Oscillating 
Mirror (ROM) was estimated to be in the range of 10 -4 — 10 -6 of the laser pulse energy for photon 
energies ranging from 20 — 40 eV, with the best results being obtained for an intermediate density 
scale-length. 
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Ultrashort XUV pulses are a promising tool for a wide 
range of applications including attosecond laser physics 
and seeding of free-electron X-ray lasers. Typically, they 
are created by the nonlinear frequency up-conversion of 
an intense femtosecond driving laser field in a gaseous 
medium. Remarkable progress has been made to the 
present date with efficiencies reaching the level of 10 -4 
at 20 nm wavelengths [U E] • Such efficiencies are not yet 
available at shorter wavelengths or for attosecond pulse 
generation and the low intensities at which harmonic con- 
version takes place in gaseous media, makes harnessing 
the high peak power in the 0.1 — 1PW regime challenging. 
High-harmonic generation at a sharp plasma-vacuum in- 
terface via the Relativistically Oscillating Mirror (ROM) 
mechanism [3] is predicted to overcome these limitations 
and result in attosecond pulses of extreme peak power 

HE]. 

While other mechanisms such as Coherent Wake Emis- 
sion (CWE) can also emit XUV harmonics [6] , the ROM 
mechanism is generally reported to dominate in the limit 
of highly relativistic intensities, where the normalized 
vector potential a 2 , = 7A 2 /(1.37 . 10 18 /im 2 W/cm 2 ) > 1. 
The efficiency of ROM harmonics is predicted to converge 
to a power law for ultra-relativistic intensities [7 , such 
that the conversion efficiency is given by 77 ~ (cj/cjo) -8 ^ 3 
up to a threshold frequency uo t ~ 7 3 , beyond which the 
spectrum decays exponentially. Here, 7 is the maximum 
value of the Lorentz-factor associated with the reflection 
point of the ROM process. While these predictions cor- 



respond well with the observations made in experiments 
using pulse durations of the order of picoseconds in terms 
of highest photon energy up to keV [HJ |9] and the slope of 
the harmonic efficiency [TO] , no absolute efficiency mea- 
surements have been reported to date. 

The plasma density scale-length plays a critical role 
in determining the response of the plasma to the inci- 
dent laser radiation. In the picosecond regime, the bal- 
ance between the laser pressure and the plasma results in 
the formation of scale-lengths and density profiles which 
are close to ideal for ROM harmonic generation in terms 
of efficiency for a broad range of laser pulse contrast. 
Achieving ultra-short (attosecond) XUV pulses requires 
lasers with 10s of femtosecond (few-cycle) duration. Un- 
der these conditions, there is insufficient time to modify 
the density scale-length significantly and hence the den- 
sity gradient and profile become critical control parame- 
ters. 

Here, we report on the first absolute measurements 
of the ROM harmonic yield. The highest yield is ob- 
served for intermediate pulse contrast, while the yield 
declines again for the highest pulse contrast, consistent 
with a plasma vacuum interface approaching step-like 
conditions. Achieving and verifying such extreme inter- 
action conditions for relativistic laser intensities is an es- 
sential step towards exploiting the potential of a wide 
range of phenomena, such as bright XUV harmonics, ra- 
diation pressure driven ion-sources and the formation of 
relativistic electron sheets [TT] . 
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FIG. 1. Experimental setup: High-contrast laser pulses were focused with an f/2 off-axis parabolic mirror to up to 3-10 19 W/cm 2 
on a fused silica or plastic coated substrate at 45° p-pol. The XUV emission was recorded with two spectrometers separately 
(in the presented data only plastic coated targets are used). The flat-field spectrometer shown in configuration 1 allows a 
measurement of the beam divergence. The XUV spectrometer system in configuration 2 has a larger collection angle and was 
calibrated regarding the incident photon flux. The black line represents the centroid beam of the laser steering into the center 
of the XUV spectrometers. 



Two experiments were performed at the 30-fs 
Titanium-Sapphire laser systems "Jeti" at the University 
of Jena and "Arcturus" at the University of Diisseldorf, 
which produced similar harmonic spectra. The laser was 
focused onto targets made of either glass or photoresist 
at an incidence angle of 45° with an FWHM intensity 
of ao = 3.5. At both laser systems the pulse contrast 
was controlled by a single plasma mirror with different 
plasma mirror targets (PMT). The plasma scale length 
L p was calculated using the hydrodynamic simulation 
code "Multi-fs" [12] based on the actual pulse profile 
measured with a 3rd order autocorrelator. Details are 
given in Ref.[13 . The highest pulse contrast was achieved 
with an anti-reflection (AR) coated PM and resulted in 
a scale length of L p < A/10, while the uncoated borosili- 
cate glass PMs produced an intermediate pulse contrast 
and L p « A/5 [14]. The harmonics' pulse energy was 
determined using an imaging XUV spectrometer with a 
8 mrad x 6 mrad acceptance angle which was calibrated 
at a synchrotron source [15] . The divergence of the har- 
monic beam was determined with an angularly resolving 
XUV spectrometer (see Fig. [l] left), which allowed the 
measured yield to be corrected for the observed angular 
divergence and thus to obtain the total photon yield. The 
measurement uncertainty of the photon yield is in the 
order of 70% (20% spectrometer calibration, 50% un- 
certainty due to divergence changes). The absolute effi- 
ciency r] = Exuv I Eo was determined by only considering 



the fraction of the laser pulse energy Eq that contributes 
to the harmonic generation process. Given the strong 
non-linearity [16 only the fraction of the laser energy 
that is focused to sufficiently high intensities (ao > 1) 
contributes to the interaction. Measurements of the in- 
tensity distribution of the focal spot with a microscope 
objective showed 20 to 40 % of the laser pulse energy 
to be concentrated within the FWHM of the focus and 
thus only this fraction is considered when comparing the 
measured efficiencies to simulations. 

Two typical harmonic spectra from photoresist targets 
are shown in Figjl] (right) using high and intermediate 
contrast settings, respectively. The angular distribution 
indicated in Figjl] (left) reveals a divergence of 18.6 mrad 
for the 21st harmonic using photoresist targets and high 
contrast. The importance of accounting for the diver- 
gence of the harmonic beam when comparing changes 
in the other parameters is highlighted by the observed 
changes in divergence under conditions where only the 
pulse contrast was varied. Changing the pulse contrast 
from the AR to the glass setting (and hence the scale- 
length from L p « A/10 to L p w A/5) changes the diver- 
gence from 18.6 mrad to 26 mrad for the 21st harmonic. 
In addition, the divergence of different ROM harmonic 
orders taken from a single measurement has an almost 
constant value. This is in excellent agreement with pre- 
vious observations [17] and fits well with the analysis that 
the divergence of the ROM harmonic beam is character- 
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istic of a beam with excellent spatial coherence and is 
determined by the curvature ('dent') of the emission sur- 
face imprinted to the target by the light pressure. Since 
the velocity of the hole-boring process depends on the 
plasma density [18] the longer scale-length deforms more 
rapidly and should therefore result in a larger divergence 
in agreement with the observations. This implies that the 
divergence should be reduced substantially in the limit 
of larger spots/or shorter pulses, which would reduce the 
curvature of the dent at the peak of the pulse respectively. 

Until now, it has generally been accepted that achiev- 
ing sufficiently steep density gradients for ROM harmon- 
ics is the major challenge and hence one would expect 
the harmonic yield to increase as the prepulse level is 
reduced. For our conditions and peak intensities, the 
strongest harmonic emission is observed for intermedi- 
ate contrast settings (glass PMTs), suggesting that the 
higher contrast setting with AR PMTs has density scale 
lengths which are even shorter than the ideal value. The 
pulse energies i^xuv (efficiencies rj) for individual har- 
monics are in the order of 3 — 24 \i J (0.1 — 1 x 10 -4 ) 
for the 17th harmonic and approximately 0.3 — 2.7 /iJ 
(0.1 — 1 x 10 -5 ) for the 21st using different contrast set- 
tings. Thus, we find for the first time that we have clear 
quantitative evidence of ROM generation in the limit of 
ultra-steep scale-lengths. While the benefit of a small, 
but finite, plasma scale-length for ROM has previously 
been highlighted by simulations [16] [19] , the experiments 
performed so-far have required the highest achievable 
pulse contrast or shortest possible scale length, respec- 
tively, in order to optimize ROM efficiency and beam 

quality [SB El- 

The influence of the plasma scale length has been 
studied both for glass substrates and photoresist tar- 
gets that have been coated onto the optically polished 
glass substrate reducing the density from 2.2g/cm 3 to 
~ 1.1 g/cm 3 (or from « 400n c to « 200n c in terms of the 
critical density n c ). The harmonic emission for these high 
target densities and respective scale lengths is compara- 
ble indicating that the enhanced harmonic emission at 
intermediate scale lengths is not very sensitive for such 
high peak densities. This means that for our parame- 
ters the harmonic emission is enhanced due to the lower 
density in the plasma gradient and not by using a lower 
maximum density. Since the reflection point of the ROM 
is located near the critical density at elongated plasma 
density ramps, the ROM process is affected by the length 
of the plasma gradient instead of the maximum plasma 
density. The observed dependence of the efficiency on the 
scale length can be understood in terms of the plasma dy- 
namics as follows. First, the denser the plasma and the 
steeper the gradient, the more the electric field in the 
skin layer is reduced. Second, the "spring constant" of 
the electron plasma becomes larger for denser and steeper 
plasmas, making the ROM harder to drive to the high 
values of 7 associated with a more efficient production of 



higher harmonic orders. 
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FIG. 2. (a) Surface field E cr i t at the critical density in 
units of the incident field E , as estimated from the equa- 
tion in Ref. [20] (black dashed line) and computed exactly 
by numerical integration for an exponential gradient (blue 
line), (b) Efficiency of SHHG above the 14th order j]rom = 
JiTwo ^( CJ ) < ^ CJ /^o for do = 3.5 at different plasma scale lengths 
from a set of ID PIC simulations. Incidence was p-polarized, 
the plasma ramp is exponential up to a maximum density of 
n e — 200n c . 

To make an analytical estimate of the field at the crit- 
ical density surface we can consider a laser interacting 
at normal incidence with the target (oblique incidence 
can be treated by switching into the frame of reference, 
in which the laser is normally incident [21]). The elec- 
trons can gain kinetic energy only through the E-field 
of the laser. This field is tangential to the surface and 
is attenuated due to the skin effect. We assume for the 
moment, that the field is non-relativistic. Evaluating the 
linear wave equation, we find the threshold condition for 
it to become relativistic. For a perfectly steep plasma 
edge, it can be calculated analytically by evaluation of 
the continuity condition at the plasma edge, yielding 
|E cr it|/|E | = 2uoq/uo p . Hence, for our laser and plasma 
parameters, the field would not be relativistic for a per- 
fect step density profile. This is reflected in Fig. [2jt>) : 
In the limit L — >• there are no relativistic harmonics. 
The skin field is however enhanced due to a finite den- 
sity ramp. As a first estimate for finite ramps we may 
consider the calculation found in Ref. [20], leading to 
|E crit |/|E | « 1A(uj L/c) 1/6 at the critical density. This 
formula becomes exact for linear and extended (L ^> A) 
gradients. For steep, exponential ramps, as are expected 
in the experiments, we find the skin field by numerical 
integration of the inhomogeneous wave equation. Results 
of this computation are shown in Fig. [5|a), along with 
the simple scaling from Ref. [20] . It can be seen that even 
small scale lengths can considerably boost the skin field 
compared to the case of step-like profiles. Already at 
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L = A/20, there is practically no attenuation at the crit- 
ical density, but the field can still grow slightly for longer 
plasma scales. We further note that the simple sixth-root 
dependence calculated for a linear gradient also yields a 
reasonable estimate for the exponential gradient, only 
slightly overestimating the field in comparison to the ex- 
act result. Fig.[2|b) shows the integrated efficiency tjrom 
of ROM harmonics for the same density gradients and a 
laser amplitude ao = 3.5. As expected from the previ- 
ous considerations and Fig{2|a), it can be seen that the 
ROM efficiency rises quickly as soon as the skin field be- 
comes relativistic. For the scale lengths L > A/10, the 
integrated efficiency remains approximately constant at 
r] « 7 x 10 -3 as expected from the ROM model. While 
the reduction in the critical field explains the drop to 
very low efficiencies at very short scale-lengths, it does 
not fully explain the efficiency scaling at intermediate 
scale-lengths. 

Since the reflection point oscillates around the immo- 
bile ion background due to the driving relativistic laser 
field, another contribution must come from the restor- 
ing force due to the quasi-static field generated by the 
plasma once the electrons are driven out of equilibrium 
by the laser field. For a given mean displacement of the 
plasma electrons, the restoring force is proportional to 
the plasma density. In the limit of a step-like density 
profile, the restoring force is directly proportional to the 
maximum plasma density, while in the limit of very long 
density scale-lengths the restoring force is determined by 
the critical density. In the intermediate case of relevance 
here, the restoring force will depend in a complex fashion 
on the density scale-length, peak density and amplitude 
of the oscillation. While it is not easily possible to express 
this dependence in a closed analytical form, it is clear 
that one would expect the effective density and hence 
the restoring force to be lower for increasingly shallow 
density gradients resulting in the dependence shown in 
Fig. [2] Generally, the denser the plasma and the steeper 
the density ramp, the harder it is for the laser to drive 
large amplitude oscillations at the plasma surface. This 
in turn leads to a smaller oscillation amplitude of the 
ROM and, consequently, to a lower 7-factor. An analysis 
of the electron spring model at ultra-relativistic intensi- 
ties can be found in Ref.[22 . In agreement with our ex- 
perimental observations a trend towards higher efficiency 
for moderately long scale-length or low peak densities is 
expected. 

Under our experimental conditions (ao = 3.5) the effi- 
ciencies are still expected to be below the relativistic limit 
regime where the r\ ~ (u/uo)~ 8 ^ 3 scaling applies. Fig. [3] 
compares our experimental results to a range of efficien- 
cies predicted by ID PIC simulations. While the efficien- 
cies are broadly compatible with the range of efficien- 
cies predicted by the simulations, they appear somewhat 
lower than predictions for the nominal density gradients 
derived from the measurement of the pulse-contrast and 




harmonic number g>/g> 

FIG. 3. Experimental efficiencies (circles) are compared to 
spectral densities from ID PIC simulations (lines) for different 
plasma scale lengths (density n e = 200n c , exponential density 
profile). The experimental efficiencies have been normalized 
to a pulse energy of 250 mJ (energy that is focused to ao > 1). 
In the ultra-relativistic limit the efficiencies converge to the 
BGP power scaling r\ w (cj/cj )" 8/3 0. 



Multi-fs modelling. While, to our knowledge, Multi-fs is 
the best suited code to calculate the hydrodynamic ex- 
pansion under such conditions, the code has not been 
validated directly by measurements of the scale-length 
under such conditions. Consequently one possible ex- 
planation for the discrepancy may be that the density 
gradients are even steeper than predicted. What is clear 
both experimentally and from simulations is that the ef- 
ficiency of the ROM process depends sensitively on the 
plasma scale length. The generation of surface waves, 
which have been found in 2D simulations, induce high 
harmonic emission at angular sidebands [23]. This may 
lead to differences between the experimental results and 
the ID simulations. Another important effect that is not 
considered in our simulations is the ion motion. In fact, 
Thaury and Quere [16] have shown that the harmonics 
efficiency in simulations with mobile ions is significantly 
reduced. 

In conclusion, we have investigated harmonic gener- 
ation in the limit of ultra-steep density gradients and 
shown first experimental evidence of the absolute yield 
reducing for very steep gradients. This demonstrates that 
relativistic interactions in the limit of ultra-steep density 
gradients can be achieved by a careful control of the laser 
parameters. Harmonic efficiency is optimized for inter- 
mediate scale- lengths. Our results suggest the generation 
of intense attosecond pulse trains with pulse energies ex- 
ceeding 10 /iJ, thus paving the way towards applications 
such as nonlinear attosecond experiments or the seeding 
of free-electron lasers with surface high-harmonic radia- 
tion. 
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